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Puff  Dispersion  Sigma  Values  from  MVP4  at  VAFB 


I.  A.  Min 
29  June  2000 


Abstract 

This  report  documents  the  calculation  of  a  values  (standard  deviations  of  material 
distribution)  for  several  series  of  SFq  puffs  released  over  Vandenberg  Air  Force  Base 
(VAFB)  in  May  of  1997  as  part  of  the  Model  Validation  Program  (MVP).  Assuming  a 
Gaussian  distribution,  estimates  of  the  a  values  are  calculated  based  on  measurements 
of  maximum  visual  extent  of  the  puffs  provided  by  R.  Abernathy  [1].  Turbulence 
fluctuation  values  derived  from  the  a  values  can  be  used  for  the  evaluation  of  REEDM’s 
turbulence  algorithm. 


1  Introduction 

Turbulent  dispersion  at  launch  sites  are  of  interest  because  it  is  the  mechanism  by  which 
the  toxic  material  in  the  launch  exhaust  clouds  left  behind  during  rocket  launches  are 
spread.  The  purpose  of  the  present  investigation  is  to  validate  the  turbulence  model  in 
the  Rocket  Exhaust  Effluent  Dispersion  Model  (REEDM),  which  is  used  operationally 
to  support  launches  in  the  Eastern  Range  (Cape  Canaveral  Air  Force  Station)  and  the 
Western  Range  (Vandenberg  Air  Force  Base).  This  validation  effort  is  one  of  the  goals  of 
the  Model  Validation  Program  (MVP). 

On  the  21st  and  the  23rd  of  May  1997,  several  series  of  SFq  puffs  were  released  from  a 
blimp  located  near  SLC-4,  at  an  altitude  of  about  1000  ttl.  The  boundary  layer  height  was 
about  400  m  on  those  days,  so  the  puff  imagery  gives  information  about  turbulence  above 
the  boundary  layer.  The  typical  vertical  extent  of  the  launch  cloud  after  stabilization  is 
between  400  to  1200  m  at  Vandenberg,  so  the  puffs  were  probing  dispersion  within  the 
envelope  of  interest. 

The  raw  imagery  data  from  three  or  four  simultaneous  infrared  cameras  was  processed 
and  compiled  by  R.  Abernathy  [1]  into  estimates  of  the  maximum  along-wind  and  cross- 
wind  puff  dimensions.  Also  provided  were  the  dimensions  perpendicular  to  each  of  the 
along  and  cross-wind  dimensions,  in  the  image.  This  information  is  used  to  derive  the 
standard  deviations  of  the  material  distribution  in  the  puff  ( ax  and  ay).  The  rates  of 
increase  of  ax  and  cry  are  used  to  estimate  the  turbulence  intensities  ( i  =  au/U  or  av/U). 

The  descriptions  of  the  experiments  and  data  acquisition  are  detailed  in  reference  [1] 
and  are  not  discussed  here.  In  section  2,  a  brief  account  of  the  methodology  used  to 
calculate  a  from  the  visual  extent  of  the  puffs  is  given.  The  a  data  is  presented  in  section 
3  in  graphical  format  while  section  4  summarizes  the  data. 
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2  Data  Analysis  Methodology 


The  quantities  measured  by  the  infrared  cameras  are  the  maximum  along  and  cross-wind 
dimensions  of  the  puffs  at  various  times  during  their  dispersion.  Also  recorded  are  the 
dimensions  perpendicular  to  each  maximum  dimension.  The  quantities  of  interest  in 
dispersion  modeling  are  the  standard  deviations  of  the  concentration  distribution  in  the 
along  and  cross-wind  directions  ( ax  and  ay).  The  vertical  concentration  distribution  ( az ) 
could  not  be  ascertained  with  sufficient  confidence  for  these  experiments  because  of  the 
camera  angles. 

Traditionally,  estimates  of  the  standard  deviations  of  the  material  distribution  in  a 
puff  are  based  [10,  8,  6,  4,  7]  on  photographic  images  of  the  line-integrated  puff  material 
(usually  smoke  particles).  Assuming  that  the  distribution  of  the  material  in  a  puff  is 
Gaussian,  the  value  of  the  integrated  (column)  density  at  the  visible  outline  is  estimated 
by  the  use  of  the  maximum  length  (or  width)  attained  by  image  over  time.  Hogstrom  [5] 
extends  this  technique  to  allow  the  observation  of  both  the  length  and  height  of  the  puff. 


The  technique  described  in  this  report  improves  upon  the  traditional  approach  in  two 
ways.  The  first  is  the  use  of  multiple  cameras  (2,  3,  or  4)  to  get  simultaneous  images  of  the 
puffs  from  different  angles.  (Frenkiel  and  Katz  [3]  also  used  two  cameras.)  Simultaneous 
pairwise  images  are  considered  and  filtered  for  acceptability  for  consideration  as  along  or 
cross-wind  view.  This  technique  allows  for  improved  determination  of  the  puff  location, 
as  well  as  puff  orientation. 


The  second  improvement  is  the  use  of  infrared  cameras  to  image  the  SF$  puffs,  which 
allows  the  use  of  the  pixel  intensity  information  to  calculate  the  actual  value  of  the  the  line- 
integrated  concentration  (column  density)  of  the  SF&,  as  formulated  by  Polak  and  Knudt- 
son  [9].  This  technique  eliminates  the  need  to  peg  the  estimation  of  the  line-integrated 
material  density  to  the  maximum  radius  achieved  by  the  cloud. 


Figures  1  and  2  show  the  distribution  of  the  column  density  or  line-of-sight  integrated 
concentration  values  for  a  puff  image  (Series  3,  Puff  1,  May  21).  The  column  density  is 
related  to  the  pixel  brightness  in  the  infrared  image,  which  was  processed  and  provided 
by  Abernathy  [2].  A  relatively  sharp  delineation  of  the  puff  boundary  can  be  seen  in  these 
figures,  imposed  by  the  noise  in  the  image  at  a  column  density  of  about  1  ppm-m.  Figures 
3  and  4  show  similar  surface  and  contour  plots  for  a  different  image  (Series  1,  PUff  1,  May 
21).  A  smoothing  algorithm  was  applied  to  these  figures  for  improved  visibility. 


Assuming  a  Gaussian  distribution  of  the  material  (in  all  three  directions),  and  assuming 
an  elliptical  shape  (so  that  when  y  =  Da/2,  z  =  0  on  the  To  =  constant  contour),  the 
equation  for  To, 


rn  = 


Q 


2iuj2a  Dperjp/  Da 


-(Da/2)*/2c\ 


(1) 


where  Q  is  the  mass  of  SF§  in  the  puff,  Da,  Dc,  and  Dperp  are  the  maximum  along-wind, 
cross-wind,  and  perperdicular  visible  extent  of  the  puff,  respectively.  To  is  the  column 
density  at  the  maximum  visible  extent  of  the  puff.  Eq  (1)  can  be  solved  using  Newton- 
Raphson  method  for  a  a-  A  solution  for  <jc  is  obtained  similarly  using 


r0  = 


Q  e-(Dc/ 

2-iralDperp/Dc 


(2) 


4 


>650 
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Figure  3:  Column  density  (ppm  —  m )  versus  cross-wind  and  vertical  dimensions  (m  x  m) 
for  May  21,  Series  1  Puff  1  image  from  camera  2  at  t=650  s 


11120650  smooth— 2 


Figure  4:  Column  density  (ppm  —  m)  versus  cross-wind  and  vertical  dimension  (m  x  m) 
for  May  21,  Series  1  Puff  1  image  from  camera  2  at  t=650  s 
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Figure  5:  Integrated  column  density  ( ppm  —  m 2)  versus  along-wind  (m)  for  May  21,  Series 
3  Puff  1  image  from  camera  3  at  t=150  s 


To  is  estimated  to  have  a  value  of  about  1  ppm-m  from  examination  of  the  processed 
infrared  images.  There  is  actually  significant  variation  in  the  exact  value  of  To,  ranging 
from  less  than  0.5  to  greater  than  2,  due  to  the  continually  varying  line-of-sight  and  relative 
position  with  respect  to  neighboring  puffs.  The  calculated  o  is  somewhat  sensitive  to  the 
value  of  To  assumed,  but  1  ppm-m  is  a  reasonable  compromise  that  bounds  the  possible 
error  to  less  than  ±30%. 

Figure  5  shows  the  result  when  the  the  column  density  values  are  integrated  along 
the  vertical  direction.  The  profile  is  compared  against  Gaussian  curves  for  two  different 
assumptions  of  the  mass  contained  in  the  puff  using  1)  Q  =  Ma  (average  puff  SFq  mass) 
and  2)  Q  =  Mj (imagery-derived  mass).  It  can  be  seen  that  the  Gaussian  assumption  is 
reasonable.  The  Gaussian  assumption  is  best  applied  for  ensemble  averages  and  is  less 
valid  for  single  realizations.  Figure  6  shows  a  similar  plot  for  a  different  image. 

Ideally,  ac  and  a  a  should  be  calculated  by  direct  integration  of  the  column  densities 
derived  from  the  pixel  intensities.  It  is  not  feasible  to  do  the  detailed  analysis  of  each  image 
required  for  the  direct  integration;  however,  it  was  done  for  several  images  to  validate  the 
approximate  Gaussian  approach.  Figures  7  and  8  show  comparisons  of  data  from  the  two 
techniques  for  Series  1  Puff  1  and  Series  2  Puff  1,  both  from  May  21.  The  figures  show 
that  there  is  reasonable  agreeement  in  the  estimate  of  a.  Specific  discrepancies  can  be 
attributed  to  several  possible  sources,  including  the  estimate  of  To  at  the  edge,  the  loss 
of  mass  outside  the  detectable  edge,  and  the  Gaussian  assumption.  Abernathy  [1]  has 
observed  that  although  there  are  significant  variations  from  puff  to  puff,  the  extent  of  the 
puff  derived  from  the  imagery  included  the  bulk  of  the  SF&  tracer. 
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Figure  6:  Integrated  column  density  {ppm  —  m2)  versus  cross-wind  (m)  for  May  21,  Series 
1  Puff  1  image  from  camera  2  at  t=650  s 
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Figure  7:  Comparison  of  cross-wind  <jq  derived  from  direct  integration  and  Gaussian 
assumption,  Series  1  Puff  1,  May  21 
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Figure  8:  Comparison  of  along- wind  <7,4  derived  from  direct  integration  and  Gaussian 
assumption,  Series  2  Puff  1,  May  21 

3  Sigma  Data 


The  estimates  of  a  a  and  ac  are  presented  graphically  in  this  section.  The  time  histories 
of  a  development  for  27  puffs  (including  along  and  cross-wind  profiles)  are  shown  for  May 
21  in  Figures  10  through  35.  For  May  23,  46  puff  0  estimates  are  shown  in  Figures  38 
through  81. 
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Figure  9:  May  21  series  1  puff  1  Crosswind  a 
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Figure  10:  May  21  series  1  puff  2  Crosswind  a 


Figure  11:  May  21  series  2  puff  1  Crosswind  a 
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Figure  12:  May  21  series  2  puff  2  Crosswind  a 
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Figure  13:  May  21  series  2  puff  3  Crosswind  a 
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Figure  14:  May  21  series  3  puff  1  Crosswind  a 
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21  s3p2  Crosswind  ac 


Figure  15:  May  21  series  3  puff  2  Crosswind  a 
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Figure  16:  May  21  series  3  puff  3  Crosswind  a 
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Figure  17:  May  21  series  4  puff  1  Crosswind  cr 


21  s4  p2  Crosswind  gc 


Figure  18:  May  21  series  4  puff  2  Crosswind  cr 
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Figure  19:  May  21  series  4  puff  3  Crosswind  a 
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Figure  20:  May  21  series  5  puff  1  Crosswind  a 


21  s5  p2  Crosswind  o( 


Figure  21:  May  21  series  5  puff  2  Crosswind  a 


21  s5  p3  Crosswind  <jc 


Figure  22:  May  21  series  5  puff  3  Crosswind  a 
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21  s2  p2  Alongwind  a 


Figure  25:  May  21  series  2  puff  2  Alongwind  a 
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Figure  26:  May  21  series  2  puff  3  Alongwind  cr 
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Figure  29:  May  21  series  3  puff  3  Alongwind  a 
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Figure  30:  May  21  series  4  puff  1  Alongwind  a 


2 1  s4  p2  Alongwind  oA 


Figure  31:  May  21  series  4  puff  2  Alongwind  a 
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Figure  32:  May  21  series  4  puff  3  Alongwind  a 
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Figure  36:  May  23  series  1  puff  1  Crosswind  a 


23  s7  p2  Crosswind  oc 


Figure  46:  May  23  series  7  puff  2  Crosswind  <7 
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Figure  47:  May  23  series  7  puff  3  Crosswind  o 
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Figure  48:  May  23  series  8  puff  1  Crosswind 
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Figure  49:  May  23  series  8  puff  2  Crosswind  a 
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Figure  50:  May  23  series  8  puff  3  Crosswind  a 
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Figure  51:  May  23  series  f  puff  1  Crosswind  a 
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Figure  52:  May  23  series  f  puff  2  Crosswind  a 
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Figure  53:  May  23  series  f  puff  3  Crosswind  a 
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Figure  56:  May  23  series  g  puff  3  Crosswind 
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Figure  57:  May  23  series  h  puff  1  Crosswind  a 


Figure  58:  May  23  series  h  puff  2  Crosswind  a 


Figure  59:  May  23  series  h  puff  3  Crosswind 
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Figure  66:  May  23  series  3  puff  1  Alongwind 
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Figure  67:  May  23  series  7  puff  1  Alongwind 


Figure  68:  May  23  series  7  puff  2  Alongwind 
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Figure  69:  May  23  series  7  puff  3  Alongwind 


Figure  72:  May  23  series  8  puff  3  Alongwind 
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Figure  73:  May  23  series  b  puff  2  Alongwind  a 


Figure  74:  May  23  series  f  puff  2  Alongwind 
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Figure  75:  May  23  series  f  puff  3  Alongwind 
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Figure  82:  Summary  of  puff  data  from  May  21 


4  Data  Summary 


Summaries  of  puff  data  are  shown  in  Figures  82  and  83  for  May  21  and  23,  respectively. 
The  derivative  da/dt  is  used  to  estimate  the  turbulent  intensity  (i  =  (Tu/U  or  i  =  av/U ), 
assuming  the  following  linear  relationship  from  Smith  and  Hay  [11] 


(3) 


for  puff  expansion.  Emperical  evidence  [8]  suggests  that  (3i  is  roughly  constant  with  a 
value  of  0.44  (assumed  here).  An  assumed  linear  relationship  between  the  turbulence 
value  and  the  rate  of  puff  a  growth  is  appropriate  if  the  time  period  of  interest  (~  500s)  of 
the  puff  growth  is  much  less  than  the  integral  time  scale  ( Ti )  of  the  atmosphere  above  the 
boundary  layer.  The  plots  indicate  that  the  linear  growth  can  be  assumed  in  many  cases. 
In  other  cases  it  is  clearly  not  appropriate  over  the  entire  regime  under  consideration. 
The  data  thus  derived  shows  turbulence  intensity  i  value  of  about  1  to  2  degrees  for  cases 
without  persistent  large  scale  shear.  These  issues  and  others,  including  comparision  with 
meteorological  data  and  application  to  REEDM  algorithm  evaluation,  will  be  discussed  in 
a  subsequent  report. 
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Figure  83:  Summary  of  puff  data  from  May  23 
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